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ON THE CAUSALITY OF PLASMA PERMITTIVITY 


ABSTRACT 


Causal solutions for effective relative complex 
permittivities of plasmas are derived. Both lossless and lossy 
plasmas are considered. 



INTRODUCTION 


The analysis of aperture antennas in the presence of plasma 
media is usually approached by assuming steady-state plane wave 
incidence on a plane plasma layer. For inhomogeneous plasma, the 
differential equations describing the fields within the plasma 
must be solved numerically in most instances. The inhomogeneity 
is expressed through the effective complex permittivity of the 
plasma media [1-3]. The permittivity is derived from the equation 
of motion of the electrons within the media for an applied 
electric field [1,3], The permittivity depends on the operating, 
collision, and plasma frequencies. 

Whenever the initial time of the input is important or a 
spectrum of frequencies is considered, one must be concerned with 
causality; i.e., the response (electric flux) of a system 
(permittivity) cannot precede the input (electric field) [4]. The 
effective complex permittivity for an inhomogeneous plasma as 
derived in the literature, which is based on steady-state plane 
wave analysis, can be easily shown to be noncausal. 

In this paper a causal solution for plasma permittivity is 
derived. The solutions are determined by using familiar 
mathematical transforms. Causality solutions for both lossless 
and lossy plasmas are considered. 
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NONCAUSAL AND CAUSAL PERMITTIVITIES 


The effective relative complex permittivity for a plasma 
media is given as [1,2] 
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where td is 
P 


the plasma frequency, v 


is the collision frequency 


(corresponds to losses) , and td is the operating frequency. For the 
permittivity given by equation (1) to be causal, its inverse 
Fourier transform must be zero for times less than zero. The 
inverse Fourier transform is given as 
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With the expression for e (td) substituted, equation (2) becomes 
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where 5(t) is the impulse function. For times greater than zero, 
equation (3) is evaluated by contour integration and becomes 
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where c > (t) represents c(t) for times greater than zero. 


Evaluating equation (3) 
yields 
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where c < (t) represents e(t) for times less than zero. The total 

function e(t) is sketched in figure 1. Since e(t) exists for times 
less than zero, the function e(u) is not causal. However, by 
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adding the quantity n — 2 5 (<d) to the right-hand side of equation 

(1), the system function can be made causal; i.e., the equation 
now representing e (id) becomes 
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The inverse Fourier transform of equation (6) readily verifies 
that the function g (cu) is indeed causal; i.e., 


c(t) = S(t) 
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which is sketched in figure 2. The effective relative complex 
permittivity for a plasma with losses of a causal system is, 
therefore, given by equation (6). 


The additional term needed for causality can be determined in 
a more rigorous way. If equation (1) is derived first in terms of 
the Laplace transform variable s, the causality solution given in 
equation (6) follows directly upon applications of the inverse 
Laplace transform and the Fourier transform. The equivalent 
expression for equation (1) in terms of the Laplace transform 
variable s is 


e(s) 


1 + 


A 

s (s+v) 


( 8 ) 


The inverse Laplace transform is written as 
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which when evaluated becomes 


c(t) = S(t) 
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Equation (10) is identical to equation (7) which has been shown to 
produce a causal result. The frequency representation given by 
equation (6) is found by taking the Fourier transform of equation 
(10) ; i.e. , 
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With equation (10) substituted, 
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The first and third integrals are readily evaluated so that 
equation (12) is written as 
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The remaining integral cannot be evaluated as it stands. However, 
rewriting equation (13) as 
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where 


sgn (t) 



t > 0 
t < 0 


equation (13) can be evaluated as 

2 

) 
v 

or equivalently, 


u~ u 2 

c( W ) = i + _E rr S((d) - j ^ 


£ 

V 


V+JO) 


2 2 
(J (d 

c(u) = 1 - _E_ + -E H 6 ((d) 
U+IS 


(d 


- J 


11 


u 2 + v 2 


(15) 
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which is identical to equation (6) . 

Thus far only permittivity with losses has been considered. 
To examine the lossless case (v = 0) , one cannot simply evaluate 
equation (17) with v = 0. However, one can substitute v = 0 in 
equation (8) yielding 
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The inverse Laplace transform of equation (18) becomes 


e(t) = 5(t) + u p *t 
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which is sketched in figure 3. The Fourier transform of 
equation (19) is 
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which cannot be evaluated in its present form. Rewriting the t 
part of the integrand as 
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where sgn(t) has been defined earlier in equation (15 ) , the 
expression for e (u) becomes 
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where the prime denotes the derivative with respect to 
u, equation (22) becomes 
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Since the Fourier transform of the product of two functions is 
proportional to the convolution of their transforms, equation (25) 
can be written as 
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where the symbol * denotes the convolution operation and ? denotes 
the Fourier transform. With the transforms substituted, equation 
(26) can be written explicitly as 
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which upon evaluation becomes 



(28) 


Equation (28) , therefore, represents the effective relative 
complex permittivity for plasmas without losses for a causal 
system. 
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CONCLUDING REMARKS 


By using familiar mathematical transforms, causal 
representations for the effective relative complex permittivities 
of plasmas have been derived. The causal permittivities for lossy 
and lossless plasmas were found to be highly singular at zero 
frequency. Permittivities which obey the principle of causality 
could be an important consideration in time-domain numerical 
solutions of Maxwell's equations. 
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t, TIME 

Fig. 1 Noncausal impulse response for plasma with losses. 
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Fig. 3 Causal impulse response for plasma without losses. 
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